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, , Abstract 

a . 

"^ ' In one-dimensional (ID) conductors with linear E-k dispersion (Dirac systems) intrabranch 

D . thermalization is favored by elastic electron-electron interaction in contrast to electron systems with 



C^ 



CJ 



a nonlinear (parabolic) dispersion. We show that under external electric fields or thermal gradients 
the carrier populations of different branches, treated as Fermi gases, have different temperatures as 



'^ ■ a consequence of self-consistent carrier- heat transport. Specifically, in the presence of elastic phonon 

O . scattering, the Wiedemann-Franz law is restricted to each branch with its specific temperature and 



is characterized by twice the Lorenz number. In addition thermoelectric power vanishes due to 



^ ■ electron-hole symmetry, which is validated by experiment. 

t^ ■ PACS numbers: 73.63.Fg, 73.23.-b, 65.80.+n 
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Within the last decades one dimensional (ID) conductors such as nanowires, nanotubes 

n 

and molecular chains have become experimentally accessible [1| . While electron populations 
at low ternperatures in ideal ID systems are predicted to behave as Tomonaga-Luttinger 
liquids [2, y] transport experiments in ID conductors have revealed various behaviors de- 
pending on the temperature range and quality of the samples. Indeed, at low temperatures 
conductance quantization \^ , and signatures of Tomonaga-Luttinger liquid 5|, Isl and Wigner 
crystallization [7| have been observed. However, the experimental realization of such systems 
is tremendously challenging and still requires further unambiguous confirmation. As temper- 
ature is increased, the features of Tomonaga-Luttinger liquids are smeared out by thermal 
broadening and carriers behave as Fermi gases. In this regime, electron transport, ranging 
from ballistic to diffusive has been successfully described by semi-classical approaches, such 



as Landauer-Biitikker formalism [8|] direct solution of the Boltzmann equation [9|, 



Monte-Carlo simulations 



10|, and 



111 ] . Paradoxically these approaches often neglect electron-electron 
(e-e) interaction as well as self-consistent heat transport regulating the energy carried by 
electrons. 

In this letter we show that in ID conductors with linear energy dispersion (Dirac system) 
energy and momentum conservation favors elastic interbranch e-e scattering, in contrast 
to ID systems with nonlinear (parabolic) dispersion[12]. As a consequence, the fermion 
populations in different branches are not in thermal equilibrium, and are characterized by 
two different temperatures, even in the lowest electric fields due to the mutual influence 
between carrier and heat transport. Our self-consistent analysis of electro-thermal transport 
of ID Dirac systems shows that the ratio between thermal and electrical conductivity is 
proportional to the branch temperature (Wiedemann- Franz law) with a factor equal to twice 
the Lorenz number. The thermoelectric power (TEP) in ID conductors vanishes because of 
electron-hole symmetry. 

The Hamiltonian describing ID Dirac systems can be written as: 

H = VFpz^y (1) 

where vp is the Fermi velocity and dy is the y component of the Pauli spin matrix, which 
gives the following E — k relationship: 

E±{k) = ±hvFk (2) 



for which the ± sign refers to two different energy branches and results in a constant density 
of states.. Here k is the wave vector along the ID z-direction. 

For linear band structure elastic e-e collisions are grouped in three classes of processes, 
i.e. intra-intra, intra-inter and inter-inter branch scattering, depending on whether the initial 
and final states remain in the same (intra) branches or change (inter) branches with colli- 
sions. Hence, we consider scattering from the initial state \ki,rii;k2,T]2) to the final state 
\k[,ri[; k'2, ri'2), where k and f] indicate the wave- vector and sign of the branch's Fermi velocity 
(?7 = +, — ), respectively (Fig. [1^). We assume for simplicity that none of these bands is 
degenerate, and that there is only one valley, but the analysis can be easily extended to 
degenerate branches and multiple valleys. We set A; = at the branches crossing (Dirac 
point), and since both momentum and energy between initial and final states are conserved, 
we get: 

ki + k2 = k[ + k'2 (3) 

E^, (ki) + Er,, (^2) = E,> (ki) + Erf^ {k2) (4) 

Because of the proportionality between E and k these two equations are linearly de- 
pendent for intra-intra branch transitions (i.e. all the electrons states belong to the same 
branch). For example, for intra-intra branch |fci, +; A;2, +) —^ |^i, +; ^2? +) scattering, multi- 
ple values of /c2 and k'2 satisfy energy and momentum conservation, given an arbitrary pair 
of values for ki and k'^ (Fig. Wp)- If any inter-branch transition takes place, Eqs. [3] and H] 
become linearly independent (as in particle systems with nonlinear E-k dispersion). Fig. [T]^ 
shows the two possible cases of inter-inter branch scattering for which the first electron tran- 
sition is |fci, — ) -^ |fc^, +). The solid (dashed) arrow to the left corresponds to the exchange 
(symmetric) scattering for which /c2 = ^i and k'2 = ki {k2 = —ki and k'2 = —k'{) which is 



totally inefficient for thermalization 12||. Intra-inter band transitions occurs if and only if 
the wave-vector of the state scattering to the different band is at the Dirac point (i.e. /c = 0) 
as shown in Fig. [T]i for |A;i,+;0,— ) -^ |/c^, +; fcg, +). Therefore, the number of collisions 
involving only intra-branch scattering processes are roughly A^^ (number of fc-states) larger 
than the number of inter-intra or inter-inter branch transitions. Since Nk is proportional 
to the number of atoms in the system {Nk ^ 1) intra-intra-branch e-e collisions are much 
more likely to occur than inter-intra or inter-inter-branch collisions. 

Because intra-intra branch scattering is more efficient than inter-intra or inter-inter 



branch scattering, carrier populations in different branches behave as independent Fermi 
gases with specific temperatures. In the presence of an electric field F an imbalance arises 
between carrier populations in different branches (different quasi-Fermi level /i^) with a 
nonzero current flow. Consequently the distribution functions in each (thermalized) branch 
reads: 

^'^^^ ^ exp[(E-/i,)ABT,] + l ^^^ 

where T^ = T^j{z) is the local electronic temperature of the branch 77. As the quasi-Fermi 

levels are far away from the band edges, and because the thermal broadening is much smaller 

than the band width, the carrier densities are independent of the respective temperatures 

due to the constant density of states. Since all the carriers share the same group velocity, 

the current is proportional to the quasi-Fermi level difference [9|], i.e. 

I = evF{n+-n_) = -^{^i+- ^iS)= g^Go^^ _ (6) 

Tin e 

where Go = e^/(7r/i) is the quantum conductance. The factor g^ accounts for the band 
degeneracy (spin degeneracy has already been considered). In Eq. [6] the current can be 
interpreted as the superposition of the electron current (/+) and hole current (/-): 

I^ = ±^{^^^-^^) = l- (7) 

where /i = (/i+ — /U_)/2 is the effective Fermi level of the system. Similarly the internal 
energy flow per unit length for each branch with respect to the effective Fermi level is: 
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^(8) 

where 6(x) is the Heaviside step function |l3l|. 

In the high (room) temperature semi-classical regime, the distribution function in each 
branch follows the stationary Boltzmann transport equation (BTE) which after making use 
of Eq. [2] reads: 

± VF [dJ±{E) + eFdEf±{E)] = dtU I ^^,,- (9) 

Here F and dtf±\^^ii are the electric field and the collision integral accounting for carrier 
scattering, respectively. We solve this equation using the method of moments, i.e. we mul- 
tiply the BTE by {E — /i)™, where m is the moment index, and integrate over energy, for 
which we obtain: 



e(F--9,/i± =± ^o',± (10) 

(^ J gc 



and 

± 5.f/± t\{f~ ia,/i,) = ^r,± (11) 

for the zeroth and first moment equations, respectively. J^™" = 2-^hv I^ooi-^ ~ t^)"^^tfri\coiidE 
is the generahzed ?7ith moment of the colhsion integral for each branch. The left-hand side of 
Eq. [TO]is the measured effective field {F' = F — dz^i±/e) |lj], and we note that in general the 
expressions for Tm^ri depend on both the local temperatures and current level (or difference 
between quasi- Fermi levels). This dependence couples the equations describing the heat 
fiow (Eq. [TTl) for each carrier population. We solve both equations for the field and the 
temperatures profiles by using the current level as a parameter. 

In the case of ID metals, the main scattering process (other than e-e interaction) is 
scattering with phonons. We compute the collision integral for this mechanism as: 



dt 



,„, = E {Rem{q) [/,'(e+) (1 - UE)) - UE) (1 - /,,(e_))] 

RaM [/v(e-)(l - UE) - UE){1 - frf{e+)]} (12) 



where e± = E ± huq and q labels the phonon wave-vector. The prefactor Rem {Rab) is 
the phonon emission (absorption) rate. Here we focus on acoustic phonon scattering {uq = 
Vs\q\), which we assume to be elastic (i.e. Vg/vp <^ 1), and for which only the longitudinal 
modes contribute to the integral. Since phonon energies are much smaller than the lattice 
temperature Tl {huq <C ksTi) and using the deformation potential approximation, we set 
Rem = Rab oc Tl [15|. Wc then define the mean free path Xac = vp/Rem- 

By subtracting and adding the 0th moments of the BTE (Eq. [TOl) for the two carrier 
branches, we obtain: 

9,(/i+-/i_) = (13) 

and 

/ = gcGoXacF' (14) 

respectively. Eq. [T^ expresses the local current conservation in the system {I (x fi^ — fi^). 
Eq. [m is nothing but Ohm's law for which the linear conductivity is inversely proportional 
to the lattice temperature {Xac oc T[^) but independent of the carrier temperatures T±. 

Using Eq.[7]we define the branch electrical conductivity as a^ = gcGoXacf^- We emphasize 
that owing to the electron- hole symmetry Eq. [TUdoes not depend on the thermal gradient. 



Therefore, the thermoelectric power (TEP) in ID (Dirac) conductors vanishes, which is in 
agreement with recent experiments on metalhc carbon nanotubes that exhibit much smaller 
TEP than that of semiconducting carbon nanotubes [l6| . 

In the presence of acoustic phonon scattering the 1st moment equation (Eq. [TTl) for each 

branch reads: 

IF' 1 
d.U± = ^ T ^ (f/+ + f/-) (15) 

Using Eqs. [8] and [131 the LHS of Eq. [15] is proportional to dzT±. This variation in the 
carrier temperature profiles is attributed to the Joule heating (first term on the RHS) and 
the inter-branch carrier scattering due to the carriers thermal imbalance (second term). 
Combining Eq. [13 for the different branches we establish the heat flow conservation: 

9,(f/+ + t/_) = IF' (16) 

which shows that all the heat is transported by the electrons and results in an inhomogeneous 
temperature profile. This is consistent with our approximation of elastic phonon scattering 
for which no energy gained by the carriers from the external field is transferred to the lattice 
(i.e. the lattice temperature remains constant). The energy production/dissipation in the 
system couples the two carrier temperatures T± (Eq. [15]). Furthermore, Eq. [T6l validates 
the assumption of the two temperature model for the electronic population in each branch. 
Indeed, if Tj^{z) = T-{z) the LHS of Eq. [T6] vanishes, which is inconsistent with a non-zero 
current. 

Substituting Eq. [8] into Eq. [16] we obtain for the heat flow: 

^acQcT^k'i.Tx. „ ^ ^^ ^^ ^acIF' , „^ 

T 3;, ^^^± = U+ + U_T ^^ (17) 

for which the pre-factor in the temperature gradient is the carrier thermal conductivity: 

f^ri — TT — '^acdc^th {i-O) 



where Gth is the quantized thermal conductance jlTI]. 
trical conductivity for each branch: 



The ratio between thermal and elec- 



K 9-772^2 



cT j-i-t Y) fJc^ 



obeys Wiedemann-Franz law with their specific carrier temperature and a proportionality 
factor that is twice the Lorenz number (for 2 and 3 dimensions) [14] . 



In the case of isothermal lattices, \ac and F are constant, and by assuming ideal contacts 
(i.e. completely absorbing and thermalizing) T-t(=pL/2) = Tq , the solution of Eq. [TSlis: 



T±{z) 



1 + 



1 + 



iT^z/L 



1 ± 2z/L / / ~ 



(20) 



for which Jtj, = [gcekBTo) / {jiK) is the current associated to the lattice temperature. The 
temperature of the carriers with positive (negative) Fermi velocity has a maximum at z = \ac 
{z = —Xac)- In the quasi-ballistic limit, Xac ^ L, the carrier temperature profiles T^ (z) 
have a linear z-dependence: 



nQB 



T^"{z) 



1 + 



3L 



JtJ V^ l) 



(21) 



To ' ■ An^Xac 

In this regime, the temperature difference T± — Tq is small even for / ~ Itq because of 
the ratio L/Xac ^ 1- In the diffusive regime {Xac *^ L), the carrier temperature profiles 
have a parabolic shape and the maximum temperature difference is of the same order as Tq 
even for small current levels. In Fig. [2] we display the temperature and heat flow profiles 
corresponding to carrier populations with positive Fermi velocity (T^{z,I) and U^{z,I)) 
for ratios L/Xac = 0.05, 1 and 4, respectively. The electron-hole symmetry in this case is 
expressed by the fact that T+(z) = T_{—z) and U+{z) = —U^{—z). 

If a lattice temperature gradient exists along the conductor Tl{z) = Tq — AT z/L), the 
carrier temperature profiles are: 



TJz) 



\ 



rr. AT\2 



L 



L + Xo 



Tn - AT- 



L 



- Tn± 



AT 



(22) 



for which we have assumed Xac{z)TL{z) = XqTq (Eq. [14]) and the symmetric boundary 
conditions Tl{-L/2) = T+{-L/2) and Tl{L/2) = r_(L/2). For AT < Tq, Eq. [22 becomes: 

AT {2z T Xq) 



T±{z)=Tq- 
and the carrier temperature difference reads 

TJz)-T(z 



2 (L + Xo) 



(23) 



AT- 



A„ 



(24) 



L + Xac 

In the case of ballistic transport [Xac/L ^ 1) the temperature difference remains T+ — 
T_ = AT along the conductor, while for diffusive transport {Xac/L ^1), the temperature 
difference between branches is negligible. 



It is important to emphasize that Eqs. [TO] and [TT] are vahd both in low and high field 
regimes. However, here we only consider elastic scattering and thereby the results are valid 
in the linear regime (i.e. low current levels and huOq <^ ksTi). At high current levels, other 
scattering processes (e.g. optical phonon scattering) have to be included in the collision 
integrals. 

In conclusion because of the effective intrabranch carrier thermalization in the high tem- 
perature regime electron populations in the different branches of ID Dirac conductors behave 
as independent Fermi gases (with their respective temperatures) out of thermal equilibrium 
as a consequence of the electro-thermal flow. In the presence of acoustic phonon scattering, 
the carrier population in each energy branch follows the Wiedemann-Franz law characterized 
by twice the Lorenz number. The TEP coefficient in ID conductors vanishes as a result of 
the electron- hole symmetry. 
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FIG. 1: Scattering diagrams due to electron-electron interaction in which energy and momen- 
tum are conserved, (a) Feynmann diagram for the elastic e-e scattering, (b) Intra-intra-branch 
scattering, (c) Inter-inter-branch scattering ; (d) Intra-inter-branch scattering 



h- 




FIG. 2: Temperature profiles (T^) and heat flow U^ for carriers with positive Fermi velocity as a 
function of the current for: A^c = 20L (a,b), Xac = L (c,d), and Xac = 0.25L (e,f). 
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